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ABSTRACT 

We present Gemini spectroscopy for 21 candidate optical counterparts to X-ray sources 
discovered in the Galactic Bulge Survey (GBS). For the majority of the 21 sources, the optical 
spectroscopy establishes that they are indeed the likely counterparts. One of the criteria we 
used for the identification was the presence of an Ha emission line. The spectra of several 
sources revealed an Ha emission line only after careful subtraction of the F or G stellar 
spectral absorption lines. In a sub-class of three of these sources the residual Ha emission line 
is broad (>400 km s“^) which suggests that it is formed in an accretion disk, whereas in other 
cases the line width is such that we currently cannot determine whether the line emission is 
formed in an active star/binary or in an accretion disk. GBS source CX377 shows this hidden- 
accretion behaviour most dramatically. The previously-identified broad Ha emission of this 
source is not present in its Gemini spectra taken ^ 1 year later. However, broad emission is 
revealed after subtracting an F6 template star spectrum. The Gemini spectra of three sources 
(CX446, CX1004, and CXB2) as well as the presence of possible eclipses in light curves of 
these sources suggest that these sources are accreting binaries viewed under a high inclination. 
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1 INTRODUCTION 

Previous surveys of faint X-ray sources have been focused on the 
Galactic Center or globular clusters. While the Galactic Center Sur¬ 
vey (e.g., Muno et al. 2003) benefits from a high source density, 
the crowding and significant extinction make the optical/infrared 
follow-up necessary for classification difficult (e.g., Mauerhan et al. 
2009). 

The Galactic Bulge Survey (GBS; Jonker et al. 2011; Jonker 
et al. 2014; Paper I & II hereafter) is a multiwavelength project that 
is designed to allow optical/infrared classification of X-ray sources 
detected in the Galactic Bulge. The GBS consists of Chandra and 
multiwavelength observations of two 6° x 1° strips centered 1.5° 
above and below the Galactic plane (see Fig. 1 of Paper I), thus 
avoiding the \b\ < 1° regions with serious crowding and extinction 
problems, while still maintaining a relatively high source density. 
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The GBS utilizes Chandra observations with an exposure of 2 ks 
for each pointing; the exposure time is chosen to maximize the rel¬ 
ative numbers of low-mass X-ray binaries (LMXBs) to cataclysmic 
variables (CVs). The completed Chandra observations have de¬ 
tected 1640 unique X-ray sources (Paper II), agreeing well with 
the estimation in Paper I which also gave a break-down of the ex¬ 
pected numbers of various kinds of objects based on source density, 
expected Chandra flux limit, etc. Among the 1640 X-ray sources, 
~ 600 are expected to be CVs, including both intermediate polars 
(IPs) and non-magnetic CVs, while the number of LMXBs is ex¬ 
pected to be ^ 250 (see Table 2 of Paper I). We also expect ~600 
chromospherically active stars or binaries, e.g., RS Canum Venati- 
corum variables (RS CVn systems; Walter et al. 1980). 

The GBS combines a large sky coverage with the good sensi¬ 
tivity to faint X-ray sources and the excellent positional accuracy 
possible with Chandra. There are two main science goals of the 
GBS (see §1 of Paper I for more details): 1) constraining the nature 
of the common-envelope phase in binary evolution by comparing 
the observed number of sources with model predictions in each 
class, e.g., CVs and LMXBs; 2) measuring the mass of the com¬ 
pact objects in X-ray binaries, e.g., eclipsing quiescent black hole 
(BH) and neutron star (NS) LMXBs, to investigate the Galactic BH 
mass distribution (e.g., Ozel et al. 2010) and to constrain the NS 
equation of state (EoS). 

Both of the science goals rely on the multi wavelength identifi¬ 
cation and classification of this large sample of faint X-ray sources. 
A variety of optical/infrared follow-up campaigns have been con¬ 
ducted. Hynes et al. (2012) identified 69 X-ray sources in the GBS 
using the Tycho-2 catalogue. These sources are coincident with or 
very close to the bright stars in that catalogue, most of which are 
likely to be the real optical counterparts to the X-ray sources. This 
sample is a mix of objects with a broad range of spectral types, in¬ 
cluding both late-type stars with coronal X-ray emission and early- 
type stars with wind X-ray emission. Many sources are foreground 
objects instead of residing in the Galactic Bulge. Britt et al. (2014) 
reported on an optical photometric survey of three quarters of the 
sky area covered by the Chandra GBS, and presented the light 
curves of variable objects consistent with the X-ray positions of 
GBS sources catalogued in Paper 1. About a quarter of the optical 
counterparts are variable, and they are expected to be a mix of IPs, 
non-magnetic CVs, LMXBs, and RS CVns. Greiss et al. (2014) 
provided likely near-infrared identification of GBS X-ray sources 
using current near-infrared sky surveys. Maccarone et al. (2012) 
found 12 candidate radio matches to the GBS X-ray sources using 
the archival NRAO VLA Sky Survey (NVSS; Condon et al. 1998). 
The majority of them appear to be background active galactic nu¬ 
clei. 

Optical/infrared spectroscopy of the detected X-ray sources is 
an essential tool to investigate their nature. Accreting binaries can 
be identified by the emission features in their optical spectra. The 
only Arm way to distinguish white dwarf (WD), NS, and BH as the 
primaries of the X-ray binaries is via measurements of the accretor 
masses, which requires high-quality optical/infrared spectroscopy. 
Britt et al. (2013) presented five accreting binaries identified in the 
GBS based on the strong emission lines in their spectra, including 
three likely IPs, one CV undergoing a dwarf nova outburst, and one 
likely quiescent LMXB (qLMXB). Torres et al. (2014) identified 
22 new accreting binaries via the Ha emission lines in their opti¬ 
cal spectra. They developed criteria of accreting binaries based on 
the equivalent width (EW) of Ha emission line (EW> 18 A), the 
breadth of the Ha emission line (EWHM > 400 km s“^), or the 
strength of He I /15876,6678 in case of narrow and weak Ha lines. 


There are also several extensive spectroscopic studies on individual 
GBS sources. Ratti et al. (2013) presented a dynamical analysis of 
GBS source CX93 in which they measured the mass of the compact 
primary and the companion star, and concluded that the source is a 
long orbital period CV. Hynes et al. (2014) identified a symbiotic 
X-ray binary with a carbon star companion in the Galactic Bulge 
based on the spectra of its optical counterpart. 

In this work, we present Gemini spectroscopy of 21 GBS 
X-ray sources with a better X-ray positional accuracy than the me¬ 
dian (due to low off-axis angles). These 21 sources are listed in 
Table [T] We will refer to these sources with their labels, i.e. CX (or 
CXB) IDs, introduced in the GBS source catalogue (Paper I and II). 
GBS sources listed in Paper I have the prefix of “CX”, while the 
remaining GBS sources, detected in the last quarter of the Chan¬ 
dra coverage published in Paper II, have the prefix of “CXB”. The 
sources in each catalogue were ranked by their Chandra counts, 
where CXI has the most counts among the CX sources. The ma¬ 
jority of the sources in this work are CXB sources, while previ¬ 
ous works were focused on CX sources. This paper is structured 
as follows. In ^we describe the Gemini observations and data re¬ 
duction. In 13 we present the analysis of the Gemini spectroscopy, 
including spectral classification and radial velocity analysis. In 0 
we give results and discuss each interesting object. Overall conclu¬ 
sions are summarized in |5] 


2 OBSERVATION & DATA REDUCTION 
2.1 Gemini Spectroscopy 

The list of 21 GBS sources of which we obtained Gemini spec¬ 
troscopy consists of six CX sources and 15 CXB sources. The six 
CX sources were proven to be interesting on the basis of earlier 
spectra and/or photometric variability. Eor example, three of the 
CX sources (CX377, CX446, and CXI004) have shown Ha emis¬ 
sion lines in their previous spectra (Torres et al. 2014) obtained by 
the Visible Multi-Object Spectrograph (VIMOS) mounted on the 
Very Large Telescope (VLT). Based on our follow-up strategy (i.e., 
prioritizing sources with higher positional accuracy and brighter 
in optical/infrared), the 15 CXB sources are chosen to have off- 
axis angles less than 5 arcmin in their Chandra observations and 
also have sufficient counts to allow for an accurate X-ray position 
(< !")• Optical/infrared brightness, colour and photometric vari¬ 
ability are also among the factors of sample selection. 

The finding charts of our sources are shown in Appendix A 
(see Eig. [T2UT3l l. The coordinates listed in Table[T]are for the candi¬ 
date optical counterparts, i.e., the objects for which we took Gem¬ 
ini spectroscopy. The astrometry was performed on images from 
VLT/VIMOS (CX377, CX446, and CX1004; Torres et al. 2014), 
Gemini (CXB 117), and Mosaic-II/Dark Energy Camera (DECam; 
other sources; see ^2.21 ). All of them have a RMS accuracy of 

< 0.2", while some sources (CX377, CX1004, and CXB 117) have 

< 0.1" positional accuracy (Britt et al. 2014; Torres et al. 2014). 

Optical spectroscopy of the 21 GBS sources was obtained with 
GMOS (Gemini Multi-Object Spectrograph; Davies et al. 1997) 
mounted on the Gemini-South Telescope in Chile. All the observa¬ 
tions were taken between 2012 Apr 20 and 2013 May 4 under pro¬ 
grammes GS-2012A-Q-44 and GS-2012A-Q-67 (see Table □ for 
an observation log). Nine objects have multi-epoch spectroscopy, 
while other objects have one epoch; each epoch has 900-3600 s 
integration time. The seeing of each spectroscopic observation was 
measured from the corresponding acquisition image (see the last 
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column of Table [T])- Typical seeing was around 0.1" (with a range 
of 0.5-1.3")- GMOS was operated in long-slit mode. We used the 
R400_G5325 grating (400 line mm“^) and a 0.75" slit. The two- 
dimensional spectra were binned by a factor of two in both spa¬ 
tial and spectral dimensions, resulting in a spatial dispersion of 
0.15"/pixel and a spectral dispersion of 1.36 A/pixel. The spec¬ 
tral resolution is estimated to be ^ 5 A FWHM for the sources 
that had filled the whole slit during the observation (i.e., the see¬ 
ing was greater than the slit width 0.75"), while spectral resolu¬ 
tion should be better than ^ 5 A FWHM for the sources with see¬ 
ing less than 0.75". All the Gemini/GMOS spectra were split into 
three equal parts in the wavelength dimension by detector gaps. 
The bluest part of the spectra was ignored in the analysis because 
of the low signal-to-noise ratio (S/N) in the spectra due to the ex¬ 
tinction towards our sources and the lack of arc lines in this part of 
the spectra. The middle part (wavelength range ^ 4800-6100 A) 
and the red part (wavelength range ^ 6200-7600 A) of the spectra 
were reduced and analyzed separately. The results presented in this 
work are mainly based on the analyses of the red part of the spectra 
because it is least affected by extinction. 

The Gemini/GMOS data were reduced using the FIGARO 
package implemented in the STARLINK software suite and the 
packages of PAMELA and MOLLY developed by T. Marsh. The two- 
dimensional spectra were bias-corrected and flat-fielded. The bias 
was corrected using the overscan areas of the detectors. We uti¬ 
lized the flat fields taken directly following each target observation 
for flat fielding. We fit the background on both sides of the spectra 
with a second-order polynomial and determined the background at 
the position of the spectra. The object spectra were optimally ex¬ 
tracted using the opt exp procedure in the PAMELA package (Marsh 
1989). The spectra were wavelength-calibrated in MOLLY using the 
helium-argon arc spectra which were taken either right after ob¬ 
serving the target or at the end of the night. The average arc spec¬ 
trum was used in case of multiple arc spectra in one night. The 
resulting root mean square (RMS) scatter on the wavelength cali¬ 
bration is <0.3 A. The wavelength calibration was examined using 
the skyline O I A6300.303. Offsets relative to the wavelength of 
this line have been corrected. For sources with seeing less than the 
slit width (0.75"), the centroiding uncertainty (i.e., if the source is 
not placed in the middle of the slit) may introduce a small wave¬ 
length shift (Bassa et al. 2006). This wavelength shift cannot be cor¬ 
rected by examining the sky-line wavelength as those fill the whole 
slit. However, it can potentially be assessed by checking the wave¬ 
lengths of diffuse interstellar bands (DIBs; Herbig 1995). Three 
such sources (CX84, CXB149, and CXB174) have strong DIBs at 
A5780 for which the line profiles do not deviate from Gaussians. 
We checked these features and found they have minor shifts relative 
to the rest-frame wavelength (65 km s“^ for CX84, —30 km s“^ for 
CXB149, and —75 km s“^ for CXB174 in heliocentric frame). We 
corrected the wavelength scale for these small shifts. Each spec¬ 
trum was normalized by dividing it by the result of fitting a 5th- 
order polynomial fit to the continuum. 

2.2 Optical Photometry from Mosaic-II & DECam 

Time-resolved optical photometry for all six CX sources was ob¬ 
tained with the Mosaic-II imager mounted on the Blanco 4-meter 
telescope at the Cerro Tololo Inter-American Observatory (CTIO) 
in 2010 July 8-15. Nineteen exposures in the Sloan /-band with 
an integration time of 120 s were taken on 45 overlapping fields 
to cover a nine square degree area, which contain all but seven of 
the X-ray sources identified in Paper I (see Britt et al. 2014 for 


the full description and results of this variability campaign). Typ¬ 
ical seeing for these imaging observations was around 1" (with a 
range of 0.8-3.0"). The data were reduced using the NOAO Mo¬ 
saic Pipeline (Shaw 2009). We performed differential photometry 
using Alard’s image subtraction ISIS (Alard & Lupton 1998; Alard 
2000) to obtain the changes in flux with respect to reference im¬ 
ages. The zero-point flux in the reference images was measured 
with either aperture photometry or DAOPHOT-II (Stetson 1987). 
The number of variable interlopers within the 95% error circles of 
the X-ray position is ^ 40 (^ 3%;see §3.1 of Britt et al. 2014 for 
details). 

The optical photometry information of four CX sources 
(CX84, CX377, CX446, CX1004) in this paper is shown in Tablegl 
while their light curves have been presented in literature (CX84 in 
Fig. 6 of Britt et al. 2014; CX377, CX446 and CX1004 in Fig. A2 
of Torres et al. 2014). For CX138, there are two blended sources 
in the Mosaic-II image at its X-ray position; neither of the possi¬ 
ble counterparts is variable. The counterpart to CXI39 is saturated 
in the Mosaic-II imaging. Among the four CX sources with light 
curves, only CX84 shows possibly periodic behaviour with a pe¬ 
riod of 4.67 days. However, it is worth noting that it is currently 
not possible to confirm the periodic nature of this variability be¬ 
cause the baseline of our monitoring was only 8 days, which is less 
than two full cycles. The light curves of the other three sources 
only have random flickering with an RMS~ 0.05-0.1 magnitude, 
although CX446 possibly experiences eclipse events. 

Optical photometry for five of the CXB sources in our sam¬ 
ple (CXB2, CXB64, CXB82, CXB99, and CXB 113) was obtained 
using the Dark Energy Camera (DECam) instrument mounted on 
the Blanco 4-meter telescope at CTIO on the nights of June 10 and 
11 of 2013. The average seeing for both nights was 1.3" with a 
range of 0.9-1.9". DECam provides a 2.2 x 2.2 square degree field 
of view combining 62 science CCDs, 8 focus CCDs, and 4 guid¬ 
ing CCDs with a scale of 0.27" per pixel. Eor all of our images, 
the SDSS r' filter was used with exposure times of either 2 x 90 s 
or 2 X 1 s for the faint and bright sources, respectively. The DE¬ 
Cam pipeline reduction provided the resampled images with cross¬ 
talk corrections, overscan, trimmed sections, bias subtraction, flat- 
fielding and saturation masksQ We then used standard IRAE tasks 
(including wcsctran, digiphot, and apphot) to extract magni¬ 
tudes and fluxes and to generate light curves through the use of dif¬ 
ferential photometry. Calibration of the target stars was achieved 
by using reference stars in the field of view that were contained in 
the Carlsberg Meridian 14 catalogue (Evans et al. 2002) and VizieR 
catalogue: 1/304. 

The photometry of these five sources are included in Table |2] 
while their light curves are shown in Eig.[T] Among the five GBS 
sources, CXB2 shows possible eclipsing/dipping events and an 
“outburst”. Although the outburst was towards the end of the night 
with higher airmass, visual inspection of the images confirms that 
the brightening is real. This outburst is possibly the reprocessed 
X-ray emission from a Type I X-ray burst, which is similar to the 
case of NS LMXB EXO 0748-676 (Hynes et al. 2006). The light 
curve shows a likely periodic modulation of P ^ 0.447 day (see 
the phase-folded light curve). CXB82 and CXB99 both appear to 
be long period variables with periods longer than our 2-night ob¬ 
serving run. Eor both sources, there is a steady increase in bright¬ 
ness by about ^0.1 and ^ 0.25 magnitude, respectively. Due to 


^ See DECam Data Handbook at http://www.noao.edu/meetings/decam/ 
media/DECam_Data_Handbook.pdf. 
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Figure 1. The DECam /-band light curves for CXB sources (CXB2, CXB64, CXB82, CXB99, and CXBl 13). CXB2 shows a periodic modulation of 0.447 
days. CXBl 13 has a sinusoidal modulation of 0.588 day in its light curve. The phase-folded light curves of these two objects are also included. 
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Figure 2. Gemini/GMOS spectra for six sources with Ha emission lines 
(CX446, CX1004, CXB2, CXB64, CXB99, and CXB113). The broad Ha 
emission of CX446, CXI004, and CXB2 likely originate in the accretion 
disk, while the narrow Ha and Hj8 emission of CXB64 and CXB113 are 
an indication of chromospherically active binaries. The Ha emission of 
CXB99 is narrow and weak. All the spectra are normalized to unity. The 
positions of the Ha and the Hj8 emission lines are labeled by the dotted 
lines. 


the lack of r' standard stars in their fields, CXB64 and CXB113 
only have approximate magnitudes without absolute calibration. 
CXB64 has a counterpart USNO Bl.O star 0578—0732346 with 
the optical magnitude of R ^ 18.7 and I = 16.9. CXB113 shows 
a sinusoidal modulation with period of P = 0.58791(12) day. The 
phase-folded light curve for this source is also shown in Fig.[T] The 
optical counterpart to CXBl 13 was also identified by OGLE (Opti¬ 
cal Gravitational Lensing Experiment; Eield: BUL_SC37, StarlD: 
9614; Udalski et al. 1992). The OGLE source has optical magni¬ 
tudes of V = 19.135(154) and/= 15.847(49). It also shows ape¬ 
riodic modulation of P 0.588 day. 


2.3 Near-Infrared Photometry 

The likely near-infrared matches to the GBS X-ray sources de¬ 
tected by Chandra (Greiss et al. 2014) were mainly obtained from 
the public variability survey VISTA Variables in the Vfa Lactea 
(VVV; Minniti et al. 2010). The VVV survey provides the most 
complete near-infrared coverage of the GBS area with consistent 
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Figure 4. Gemini/GMOS spectra for the remaining 6 sources presented in 
this paper. All the spectra are normalized to unity. The positions of the Ha 
and the Hj8 emission lines are labeled by the dotted lines. 


depth. Complementary coverage from the Two-Micron All Sky 
Survey (2MASS; Skrutskie et al. 2006) and UKIDSS Galactic 
Plane Survey (GPS; Lucas et al. 2008) were also utilized for bright 
{Ks < 12.5) and faint {Kg > 16) sources, respectively. Greiss et al. 
(2014) developed a method to estimate the likelihood (by calculat¬ 
ing the false alarm probability) of a near-infrared source to be the 
counterpart of the X-ray source. However, for the optical sources 
where we have obtained Gemini/GMOS spectra, we are able to 
search for the near-infrared counterparts that matches our optical 
sources by comparing the optical and near-infrared images. We 
set a 0.2" astrometric error circle, and searched for any matches 
within that error circle between the optical images and VVV im¬ 
ages. Then we visually examined both images and selected the true 
near-infrared counterparts. A visual inspection is crucial given that 
the GBS fields are crowded. 

Table m lists the near-infrared counterparts to 13 of the 21 
optical sources with Gemini/GMOS spectroscopy. The other eight 
sources do not have reliable near-infrared counterparts. Three of 
them (CX446, CXB26, and CXB137) are too faint in the VVV 
frames and there are no near-infrared sources at the optical position. 
For the other five sources (CX1004, CXB64, CXBl 17, CXB189, 
and CXB201), there could be near-infrared matches but they are 
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Figure 3. Gemini/GMOS spectra for nine sources with Ha absorption features (CX138, CX377, CXB26, CXB73, CXB117, CXB149, CXB174, CXB189, 
and CXB201), many of which also have Hj8 absorption. All the spectra are normalized to unity. The positions of the Ha and the Hj8 emission lines are labeled 
by the dotted lines. 


blended with nearby sources, for which the current VVV data re¬ 
lease does not provide photometry. 


The infrared colours can be used to estimate the distance of the 
source by calculating the needed extinction to match the source in¬ 
frared colours to those of a standard star with the same spectral type 
(e.g., Ratti et al. 2013). Some of the GBS sources have an infrared 
excess, i.e., their infrared colours are redder than expected based 
on optical colours (see §4 of Hynes et al. 2012), which is possibly 
for cases where the optical classification may represent the hotter 
component of a binary (or of a blend, if another star is serendipi- 
tously along the same line of sight) while the infrared colours are 
from the cooler component. Accreting systems (CVs and LMXBs) 
and Be stars with a circumstellar disk may also have colours that do 
not match those of single stars. Be stars could also have significant 
colour variations due to the formation/dissolution of the circum¬ 
stellar disks. 


3 DATA ANALYSIS 

The Gemini/GMOS spectra of the 21 GBS sources are shown in 
Fig-EHU Three sources (CX446, CX1004, and CXB2) show appar¬ 
ent broad Ha emission lines, while another three (CXB64, CXB99, 
and CXB113) have narrow Ha emission (see Fig.|2]). We fit Gaus¬ 
sian profiles to these broad Ha emission lines to measure their 
width and velocity separations (if double-peaked) using the mgf it 
procedure in the MOLLY package. The equivalent widths (EWs) 
of these Ha lines were measured using the light procedure in 
MOLLY. Nine other objects have Ha absorption features (Fig. [3]). 
The stellar absorption features shown in the Gemini/GMOS spectra 
can be utilized to perform spectral classification, and radial veloc¬ 
ity analysis. The spectra of the remaining six sources are shown in 
Fig. SI they appear at first sight to have neither Ha emission fea¬ 
tures nor Ha absorption features. 

Some of the objects in our sample have strong DIBs in their 
spectra. We estimated the reddening for these sources via the equiv¬ 
alent width (EW) of the DIB at A5780 with the calibration in Ta¬ 
ble 3 of Herbig (1993). We also compared the measured reddening 
to the Bulge reddening along the line of sight provided by Gonza- 
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lez et al. (2011,2012), which utilized the Red Clump stars in the 
Bulgell 


3.1 Spectral Classification 

We utilize the optimal subtraetion technique following previous 
works (e.g.. Marsh et al. 1994; Ratti et al. 2013) to classify the 
spectra. A set of template star spectra was chosen from the library 
of the Ultraviolet and Visual Echelle Spectrograph Paranal Obser¬ 
vatory Projeet (UVES POP; Jehin et al. 2005), eovering spectral 
types from AO to M6 with luminosity class of V (see Table ID). The 
UVES POP template spectra provide coverage over the wavelength 
range of 3000-10000 A with a spectral resolution of ^ 80,000. The 
templates were re-sampled and Gaussian-smoothed to match the 
spectral resolution of the object spectra. The object spectra were 
Doppler-shifted into the same rest frame and averaged. Each stel¬ 
lar template was optimally subtracted from the object spectrum, 
while a test was performed on the residuals. All the emission 
lines, DIBs, and telluric features (e.g., Kuruez 2006; Wallace et al. 
2011) were masked during the procedure. The resulting x^ values 
for each template were eompared with each other. The template 
with the minimum x^ value provides our best estimate of the spee- 
tral classification of each GBS source in our sample. These spectral 
classification procedures were first performed on the red part of 
the Gemini spectra. We verified our results by performing the same 
procedures to the middle part of the Gemini spectra and the results 
are consistent with each other. The results of spectral classifications 
are listed in Table |5] Note we are not controlling luminosity class 
in the spectral classification. Nine of the 21 GBS sources in our 
sample were spectrally classified. The uncertainty of spectral clas¬ 
sification is estimated to be one or two speetral sub-classes. The 
other three sources listed in Table [5] are classified as early-mid M- 
type based on their spectral features. 

We have also tried to measure the rotational broadening of 
the absorption features. However, our limited speetral resolution of 
^ 200 km s“^ precluded the determination of rotational velocities. 

3.2 Radial Velocity Analysis 

The radial veloeities of the optieal eounterpart of each GBS source 
are measured by cross-correlating the object spectra using the xcor 
procedure in MOLLY. Eor each source, all the object spectra and the 
template spectrum were rebinned to the same veloeity dispersion 
with the vbin procedure. We also masked all the emission lines, 
interstellar/telluric features during the cross-correlation procedure. 
We performed two sets of the cross-correlation analysis. The first 
set uses the first souree spectrum as the cross-correlation template, 
i.e., we derive the radial velocity (RVl) relative to the first source 
speetrum. For the GBS sources that were spectrally classified using 
the procedure in dl we also performed a second set of cross¬ 
correlation analyses, taking the UVES POP standard star with the 
best-fit speetral type as the template. Thus the seeond set of radial 
velocity (RV2) is relative to the template star. All spectra had been 
shifted to a heliocentric frame. 

Radial velocity values were derived by fitting a Gaussian pro¬ 
file to the eross-eorrelation function. The results are listed in Ta¬ 
ble [6] Seven sources have non-zero radial velocities relative to the 

^ See the Bulge Extinction And Metalicity (BEAM) calculator at 
http://niill.astro.puc.cl/BEAM/calculator.php 



Figure 5. The Ha region of both epochs of Gemini/GMOS spectra for 
CX446, acquired on 2012 Jun 22 (blue line) and 2013 May 4 (red line), 
respectively, showing evidence for strong variability in the EW of the Ha 
between the two epochs. The VLT/VIMOS spectrum of CX446 is also over¬ 
laid (black line). 

standard star template. Three sources (CX84, CX138, CX139) have 
radial velocity variations of ^ 100 km s“^ on timescales of days. 


4 RESULTS & DISCUSSION 

4.1 Objects with Ha Emission Lines 

There are six sources with apparent Ha emission in their Gem¬ 
ini/GMOS speetra (Fig.|2]). Three of them (CX446, CX1004, and 
CXB2) have broad Ha emission lines (FWHM>800 km s“^), 
while those of the other three sourees (CXB64, CXB99, and 
CXB113) are narrow (FWHM<200 km s“^). The strong, broad 
Ha emission likely originates from an accretion disk, indicating 
CX446, CXI004, and CXB2 are likely acereting binaries ( ^4.1.11 ). 

CXB64 and CXB113 are possibly chromospherically active 
stars or binaries beeause of their narrow and weak Ha emission 
lines (Torres et al. 2014; also see Fig.|2]). Th^ H — K eolour (0.300± 
0.033; see Table [3]) of CXB113 agrees well with that of a M4V- 
M5V star (H — K ^ 0.30; see Table 5 of Pecaut & Mamajek 2013), 
while its J — H colour (0.509 ± 0.030) appears slightly bluer than 
that of a M4V-M5V type star (J — H^ 0.57). 

CXB99 has a weak, narrow Ha emission. After the optimal 
subtraction with the best-fit K2V template, the Ha emission ap¬ 
pears to be stronger, indieating that it partially fills in the Ha ab¬ 
sorption line in the spectrum. 

4.1.1 Potential Quiescent Accreting Binaries: CX446, CX1004, 
and CXB2 

Both epochs of the Gemini/GMOS spectra of CX446 show broad 
Ha emission (see Fig. |3. The strength of the Ha line is weaker 
in the spectra taken on 2013 May 4 than that on 2012 Jun 22. The 
line width (1250 ± 50 km s“^) is smaller than that of the Ha emis¬ 
sion lines in the VLT/VIMOS spectra of CX446 (2200±50 km s”!; 
Torres et al. 2014), while the line EW is bigger. No absorption lines 
from the companion star are visible in the Gemini speetra. We also 
see no evidenee of He I A 6678. The light curve of CX446 does 
show a possible eclipse event with a depth of 0.4 magnitude (see 
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Figure 6. The Ha region of the four epochs of Gemini/GMOS spectra for 
CXI 004 (from top to bottom in chronological order; see Table [I). All four 
spectra show a broad, double-peaked Ha prohle which varied in strength 
between the observations. 
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Figure 7. The Ha region of the four epochs of Gemini/GMOS spectra for 
CXB2 (from top to bottom in chronological order; see Table[T). The spectra 
of the first and third epochs show broad, weak Ha emission, whereas it is 
absent in the bottom spectrum. Both the emission and absorption features 
have substantial negative velocity offset which may indicate the high space 
velocity of this object. 
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Fig. A2 of Torres et al. 2014; on HJD= 2455387.82). No signifi¬ 
cant periodicity is found in the light curve. CX446 is a candidate 
eclipsing CV or qLMXB. 

The Gemini/GMOS spectra of CXI004 show double-peaked 
Ha emission in all four epochs (see Fig. [6]). We measure the Ha 
emission line properties using the spectra of the last epoch because 
of its best S/N. The line width is FWHM = 2500 ± 100 km s“^, 
which is broader than that in the VLT/VIMOS spectra of CXI004 
(FWHM = 2100 ±20 km s“^ see Torres et al. 2014). The line 
strength (EW= 38.0 ± 0.6 A) has also slightly increased com¬ 
pared to that in the VLT/VIMOS spectra (EW= 32.9 ± 0.4 A). 
The velocity separation between the red and blue peaks is Av = 
1160 ±30 km s“^, which is consistent with the result in Tor¬ 
res et al. (2014). The centroid of the Ha line does not have sig¬ 
nificant radial velocity (velocity offset —15 ±20 km s“^), while 
that of the VLT/VIMOS spectra shows substantial radial velocity 


(—170±20 km s“^; Torres et al. 2014). One possible scenario for 
this is that the accretion disk of CXI004 is processing. The double- 
peaked profile of the Ha emission line is asymmetric in three of 
the four epochs; the relative strength of the blue peak and the red 
peak varies with time (see Eig. El. This behaviour can be explained 
by the presence of an S-wave originating in a hot-spot or the donor 
star (see §IV of Johnston et al. 1989). 

The prominent molecular absorptions in the wavelength range 
of 6300-7300 A are signatures of an M-type companion star. How¬ 
ever, the detailed rotational broadening analysis and spectral classi¬ 
fication are not feasible for the Gemini/GMOS spectra of CXI004. 
As discussed in Torres et al. (2014), the line width and the double¬ 
peak velocity separation are consistent with that of an eclipsing qui¬ 
escent CV or qLMXB. However, possibly owing to the faintness of 
the source, the light curve of CXI004 (see Eig. A2 of Torres et al. 
2014) does not show significant periodic behaviour. Torres et al. 
(2014) suggested CXI004 to be a nearby source due to the lack of 
diffuse interstellar bands. 

The four epochs of Gemini/GMOS spectra of CXB2 are 
shown in Eig. |7] The spectra of the first and third epochs appear to 
have weak, broad Ha emission. Both of the lines have large veloc¬ 
ity offsets (which are also visible in the figure) with Vpeak = —260± 
60 km s“^ and —200± 30 km s“\ respectively. The widths of the 
two lines are EWHM = 825 ± 175 km s“^ and 725 ± 100 km s“^, 
respectively. The spectrum of the fourth epoch shows a narrow Ha 
absorption feature (EWHM = 180 ±40 km s“^), which also has 
a significant velocity offset (vpeak = —175 ± 15 km s“^). These 
large negative offsets could originate in a high space velocity of 
this source, which may imply a NS or a BH as the primary star. 
CXB2 also shows the Ha emission variation as we have seen in 
CX377. No photospheric line from the companion star is detected 
from the spectra of CXB2. 

CXB2 is brighter in X-rays than most of the GBS sources. 
It has 147 Chandra ACIS-I counts in 0.3-8.0 keV band in 
a 2 ks exposure. CXB2 was serendipitously detected by an 
archival Suzaku observation with 53 ks exposure time (ObsID: 
507031010). Eull details of the X-ray analysis will be presented 
elsewhere, but we briefly summarize relevant results here. The 
Suzaku XISO and XISl data were reduced and analysed using 
standard FTOOLS. A good spectral fit (;^^/v=1.0) was found for 
CXB2 with an absorbed power-law model (photon index L = 
0.96 ±0.11, Vh = L5+}-4 X 10^1 cm-2). The X-ray flux is 1.5 x 
10“^^ erg cm“^ s“^, which is consistent with that from our Chan¬ 
dra observation (1.3 x 10“^^ erg cm“^ s“^ assuming the same 
spectral model). This source is probably associated with the AACA 
source AX J1754.0—2929 which is 7.2" away (see catalogue in 
Sakano et al. 2002). The X-ray flux in the ASCA observation is 
1.4 X 10“^^ cm“^ s“^ in 0.7-10 keV band, which is similar to the 
X-ray flux in the Chandra observation. This source was, however, 
not detected in ROSAT and is classified as a transient (Paper II). 
The Nil value is smaller than that predicted to the Galactic Bulge 
(~ 10^^ cm“^), hence we infer that CXB2 probably lies in the 1- 
4 kpc distance range. Its near-infrared colours also indicate low 
reddening. The X-ray luminosity would be 1-40 x 10^^ erg s“^ (as¬ 
suming the above distance range), which is consistent with either a 
high-inclination CV or a qLMXB. 


4.2 Objects with Ha in Absorption 

Nine sources (CX138, CX377, CXB26, CXB73, CXB117, 
CXB149, CXB174, CXB189, and CXB201) in our sample only 
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Figure 8. The residual spectra for the nine spectrally classified sources (CX84, CX138, CX139, CX377, CXB73, CXB82, CXB99, CXB149, and CXB174) 
after optimally subtracting the standard star template. The dotted line in each panel labels the position of Ha. All the residual spectra show Ha emission 
features, except CXB149 and CXB174. For all panels we plot the observed spectrum on top, the best-fit template star spectrum in the middle and the residual 
spectrum, obtained after subtracting the best-fit template from the observed spectrum, at the bottom. The offset between the three spectra per panel is chosen 
so that the spectral features are discernible. The interstellar absorption feature at ~ 6280 A is masked out when optimal subtraction was performed. 


show Ha absorption lines in their Gemini/GMOS spectra (Fig.O. 
We are able to obtain a spectral classification for five of them 
(CX138, CX377, CXB73, CXB149, and CXB174). The residual 
spectra of these sources after optimally subtracting their best-fit 
stellar templates are shown in Fig. [8] Three of them (CXI38, 
CX377, and CXB73) have unambiguous Ha emission features, 
which partially fill in the stellar Ha absorption features in our 
Gemini/GMOS spectra. We measured the strength (EW) and width 
(FWHM) of the Ha emission shown in these residual spectra (see 
Table [7]). The broad Ha emission (FWHM>400 km s“^; see Ta- 
ble|7]and Fig. [8]) is most likely due to the accretion disk. Therefore, 
for the first time, we may have discovered a hidden population 
of accreting binaries where the accretion disk contribution to the 
optical light is small, and the Ha emission lines become apparent 
only after subtracting the stellar contribution to the spectra. In what 
follows, we provide further details on the spectroscopic and/or 
photometric properties of several individual objects. 

(a) CXI 38 has a broad Ha emission line (FWHM= 
1350 ± 75 km s“^) in its residual spectrum. It also has a ra¬ 


dial velocity variation of ~ 100 km s“^ between its Gemini 
spectra. Therefore, CXI38 is a likely candidate of the hidden 
accreting binaries. The reddening measured via the strength of 
DIB A5780 (EW = 1971±131 mA) is Ak = 1.29 ±0.08, which 
also puts the object in the Galactic Bulge, if not farther. Eor the 
Galactic Bulge distance (^ 8 kpc), the absolute V-band magnitude 
of CXB138 would be My = —2, suggesting a giant companion star. 

(b) CX311 is also a candidate hidden accreting binary. The Gem¬ 
ini/GMOS spectra of CX377 do not show any notable emission 
features (see Eig. [3] and Eig. [9]). However, the optical spectra of 
CX377 taken ~ 1 year before (2011 May 28 and 2011 July 23) by 
VLT/VIMOS show strong broad Ha emission with an asymmetric 
double-peaked profile (Torres et al. 2014; see Eig. |9}. The intrinsic 
width of the Ha emission is ^ 1200 km s“^, while the double¬ 
peak separation is ~ 700 km s“^ This apparent dramatic change 
of the spectrum is not caused by incorrect targeting. The finding 
charts of CX377 for our Gemini/GMOS spectroscopy and for the 
VLT/VIMOS spectroscopy (Fig.[T0l also see Eig. 9 in Torres et al. 
2014) show the exact same targets. The slit was placed during our 
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Figure 9. The comparison of the Ha region of CX377 optical spectra, ac¬ 
quired by VLTA^IMOS (upper line) and Gemini/GMOS (lower line), re¬ 
spectively. The broad Ha emission line in the VLTA^IMOS spectrum is 
absent in the Gemini/GMOS spectrum, indicating strong spectral variabil¬ 
ity even though the source remains in quiescence. 

Gemini/GMOS observations (see Fig. Gil to contain both the target 
and the adjacent object so that we are able to distinguish their spec¬ 
tra. The seeing during the Gemini/GMOS observations (~ 0.6") 
was better than that during the VLT/VIMOS observations (^ 1"). 
Thus, the contamination of the nearby bright sources is less of a 
problem for the Gemini/GMOS spectra. Therefore, the Ha emis¬ 
sion feature of CX377 indeed has significant variation within one 
year. Although the Gemini/GMOS spectra of CX377 only show 
Ha absorption lines, its residual spectrum after optimally subtract¬ 
ing the F6V standard star template shows broad Ha emission (see 
Fig.U. The Ha emission has not completely disappeared, but it 
is significantly weaker during our Gemini/GMOS observation than 
during the VLT/VIMOS observations of Torres et al. (2014). It has 
also become narrower (FWHM = 660 ±30 km s“^ in the Gem¬ 
ini/GMOS residual spectrum). This rare behaviour of Ha emis¬ 
sion line variations indicates significant accretion disk variations 
for CX377. The soft X-ray transient GRO J1655—40 has shown 
similar strong variability of its Ha emission line (e.g., Soria et al. 
2000). GRO J1655—40 has a period of 2.62 days; the spectral type 
of its companion star had been classified as F3-F6 (Orosz & Bai- 
lyn 1997), similar to that of CX377 (see below). However, the Ha 
emission line variability of GRO J1655—40 occurred when this 
object was going through an outburst cycle, while CX377 has re¬ 
mained in quiescence. It is possible that CX377 experienced a faint 
X-ray outburst ending before June 1, 2012 when our Gemini spec¬ 
tra were taken, as long as the X-ray fiux remained below the detec¬ 
tion threshold of all sky monitors such as Monitor of All-sky X-ray 
Image (MAXI; Matsuoka et al. 2009). Another possibility is that 
CX377 has variable accretion rate during quiescence. 

CX377 is spectrally classified as F6 type. The relative radial 
velocities between the Gemini spectra of CX377 (see Table [6]) are 
consistent with zero or less than 10 km s“^ which is not surpris¬ 
ing given that the three object spectra were only 15 minutes apart, 
and the orbital period of CX377 is expected to be much longer than 
15 minutes for an F6 main sequence star filling its Roche lobe. 
Then we cross-correlated the object spectra to the F6V star tem¬ 
plate (HD 16673). The radial velocities relative to the star template 
range from —50 to —20 km s“^ The radial velocity of the star tem¬ 
plate itself is consistent with zero (—4 ± 5 km s“^). The light curve 
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Figure 10. The finding charts of CX377 for our Gemini/GMOS spec¬ 
troscopy (left panel) and the VLT/VIMOS spectroscopy in Torres et al. 
(2014; right panel). The two short thick bars in each panel indicate the po¬ 
sition of CX377, while the thin lines show the slit position. The sky area 
shown in both panels is 15" x 15". North is up and east is left in both pan¬ 
els. 


of CX377 obtained in July 2010 only shows non-periodic dickering 
with a RMS scatter of 0.06 magnitude (see Fig. A2 of Torres et al. 
2014). The strength of DIB A5780 (EW = 1093 ± 86 mA) in the 
spectra of CX377 corresponds to Ak = 0.73 ±0.06. Extinction map 
from Gonzalez et al. (2012) gives Ak = 0.72 ±0.15 for CX377, 
which is consistent with the extinction derived above. Therefore, it 
is likely that CX377 resides in the Galactic Bulge. 

The neutral hydrogen column density, Nh = (1.44 ±0.14) x 
10^^ cm“^, was estimated from the extinction following the rela¬ 
tion in Giiver & Ozel (2009). CX377 has 7 counts in the ACTS-I 
0.3-8 keV band in the 2 ks exposure. Assuming an absorbed 
power-law spectrum with photon index of 1.6, the unabsorbed 
X-ray flux of CX377 is ^ I x 10“^^ erg cm“^ s“^. At a distance 
of ^ 8 kpc, the X-ray luminosity would be Lx = 8 x 10^^ erg s“^. 
With the estimated extinction and the apparent /-band magnitude 
(see Table 121), we calculate the absolute magnitude of M^' = —1.3 
and the X-ray to optical flux ratio fx/fr' ^ 0.005; both values are 
consistent with those for a giant companion star. However, the disk 
luminosity at the time of the VIMOS observation should have been 
larger than that of the E6 giant (M^/ = —1.3), implying that the 
source might have experienced a faint X-ray outburst (cf. Wijnands 
& Degenaar 2013). 

(c) CXB149 and CXB174 were spectrally classified as G6 and 
E4, respectively. Their residual spectra after optimal subtraction 
do not show unambiguous Ha emission or absorption features. 
The Gemini/GMOS spectra could possibly be from interlopers 
instead of the true optical counterparts of the X-ray sources. The 
Ha absorption features of these two sources are significantly 
blueshifted from the laboratory wavelength (—109 ±7 km s“^ 
for CXB149 and -179 ± 8 km s"! for CXB174). This shift 
is not introduced by the wavelength calibration process since 
offsets to the sky lines were small and they have been corrected. 
It cannot be explained by the wavelength shifts caused by the 
centroiding uncertainty within the slit (which are —30 km s“^ 
and —75 km s“^ respectively; see HD. These velocities are 
consistent with that for stars residing in the Galactic Bulge (e.g., 
Zoccali et al. 2014). Regarding the reddening measurements, 
CXB149 has EW(DIB A5780) = 1099±28 mA; the corresponding 
reddening is E(B-V) = 2.11±0.05 and Ak = 0.73 ±0.02, which 
put this object in the Galactic Bulge or farther (Ak = 0.30 ±0.07 
from the map of Gonzalez et al. 2012). Eor CXB174, EW(DIB 
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Figure 11. The residual spectra of CX84 after subtracting the best-fit star 
template (G9), and other templates with spectra types (F6, G6, K2) brack¬ 
eting that of the best-fit template. The Ha emission feature exists in all four 
residual spectra. 

A5780) = 616±34 mA; the corresponding reddening is E(B-V) 
= L21±0.06 and = 0.42 ±0.02, which is consistent with the 
Bulge reddening (Ak = 0.35 ±0.08) from Gonzalez et al. (2012). 
The absolute magnitudes for CXB149 and CXB174 residing in the 
Bulge would be My = 0 and My = —0.7, respectively, which are 
consistent with those of the giant stars with their spectral types. 

(d) CXB26 and CXB189 were not spectrally classified using our 
method. CXB26 is considered to be likely associated with an 
OGLE source 1" away from the Chandra position (OGLE BUL- 
SC3 6033) which is identified as a CV (see Paper II). However, 
the Gemini/GMOS spectra of CXB26 do not show apparent fea¬ 
tures of an accretion disk. It is possibly a hidden accreting binary, 
which could be verified after it is spectrally classified with future 
spectroscopy. The physical nature of CXB189 remains unclear. 


4.3 Objects without Apparent Ha Emission/Absorption 
Features 

The remaining six sources show neither clear Ha emission nor ab¬ 
sorption features in their Gemini/GMOS spectra (see Eig.|4l). Three 
of them (CX84, CXI39, CXB82) were spectrally classified (see 
Table O. The residual spectra of these sources after optimally sub¬ 
tracting their best-fit stellar templates are shown in Eig. [8] Similar 
to CXI38, CX377, and CXB73, they all have Ha emission fea¬ 
tures, which fill in the stellar Ha absorption features in our Gem¬ 
ini/GMOS spectra. To test the robustness of our results, we also 
subtracted other non-best-fit stellar template from the object spectra 
to see how the residual Ha emission would vary. Eig.fTTlshows one 
example of this test. The Ha emission feature remains in the resid¬ 
ual spectra of CX84 after subtracting several stellar templates with 
spectral types bracketing that of the best-fit template, which sup¬ 
ports the detection of a residual accretion disk. Note that the choice 
of templates with non-best-fit spectral types could affect the profile 
of the residual Ha emission. Eor example, the mis-matched tem¬ 


plates generate broader wings of the Ha emission which may result 
in a higher EWHM value (e.g., EWHM increases from 225 km s“^ 
to 320 km s“^ if using an E6 template spectrum). 

It is worth noting that even in the case that the residual Ha 
emission features are real, there could be alternative explanations. 
Eor example, it is possible that the observed optical light is dom¬ 
inated by a foreground, physically unrelated star (instead of the 
companion star in the above scenario) so that we can only see the 
accretion after subtracting the foreground stellar emission. How¬ 
ever, based on the analysis in Britt et al. (2014), we expect ~ 0.7 
interloper in our sample of 21 GBS sources (also see ^2.21 ). In 
^4.21 we already discussed two possible interlopers (for CXB149 
and CXB174) for which the residual spectra do not have clear Ha 
emission. Therefore, it is not likely that the six optical sources that 
only have Ha emission in their residual spectra are all chance 
interlopers. Another possibility is that these sources are active 
stars/binaries, where the Ha emission fills in the stellar absorption 
features. Sources with narrower Ha emission (e.g., <200 km s“^) 
in their residual spectra could fit into this scenario. However, it is 
difficult to explain the broad residual Ha emission lines of CXI38 
and CX377. The periodic modulation in the light curves of CX84, 
CXB82, and CXB99 could originate from either binary motion or 
stellar activity (hotspots). Alternative explanations for the narrow¬ 
ness of their Ha emission lines are that the binaries are relatively 
face-on, or that the disks are highly recessed. 

CX84 is spectrally classified as G9. The reddening estimated 
from the strength of interstellar band A5780 (EW = 2499±218 mA) 
is Ak = 1.62 ±0.14, which means the object is in the Bulge or 
farther. Assuming CX84 resides in the Galactic Bulge, the absolute 
V-band magnitude would be My = — 1, which indicates a G9 
giant star. The X-ray luminosity is ~ 2.6 x 10^^ erg s“^ which 
appears higher than that for G9 giants (see Eig. 2 of Giidel 2004). 
However, it is typical of long orbital period qLMXBs (e.g., Jonker 
et al. 2012). The optical light curve of CX84 indeed shows a 
modulation with a tentative period of 4.67 days. If it is a qLMXB, 
the orbital period should be twice that value (p:^ 9.3 days) because 
of the ellipsoidal modulation. In summary, the X-ray luminosity 
of CX84 suggests the qLMXB nature, while the narrow residual 
Ha emission line favors the chromospherically active star/binary 
scenario. Similarly, the reddening of CXB82 derived from its near- 
infrared colour {Ak ~ 0.32) is also consistent with the reddening 
expected if it resides in the Galactic Bulge. The absolute V-band 
magnitude {My = —0.3) suggests a G9 giant star, while its X-ray 
luminosity (^ 1.1 x 10^^ erg s“^) indicates that it is more likely 
to be a qLMXB. The Ha emission line in the residual spectrum 
of CXB82 is also narrow. Based on current data, we are not able 
to determine whether CX84 and CXB82 are chromospherically 
active stars/binaries or qLMXBs. 

Among the seven potential accreting binaries identified in 
our sample with a spectral classification, one source (CX377) has 
an E-type companion star, another four (CX84, CX138, CXB73, 
CXB82) have G-type companion stars, while the other two (CXI39 
and CXB99) have K-type companion stars. Given the mass of E/G- 
type companion stars, the mass ratios {q = M 2 /M\ here) for the 
sources with F/G-type companion stars should be higher than those 
of the more typical X-ray binaries with later spectral types (K, M). 
In contrast to LMXBs, very few of these “intermediate” mass X-ray 
binaries (IMXBs) have been found. Existing examples include Cyg 
X-2 (e.g., Casares et al. 2010) and Her X-1 (e.g., Reynolds et al. 
1997). Our Gemini/GMOS spectroscopy indicates that the GBS 
may be able to provide an example of this type of X-ray binary. 
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5 CONCLUSIONS 

The GBS is a multi wavelength survey project where one of the 
goals is to identify X-ray binaries in the Galactic Bulge area. In 
this work, we present optical and infrared photometry, and Gem- 
ini/GMOS spectroscopy of 21 GBS X-ray sources detected in the 
GBS. One prime goal of the GBS is to identify eclipsing qLMXBs. 
CX446, CXI004 and CXB2 are promising candidates to be eclips¬ 
ing qLMXBs. CX446 has broad Ha emission, and its optical light 
curve shows likely eclipsing events. CXI004 has broad, double- 
peaked Ha emission. The light curve of CXB2 also contains eclips¬ 
ing events suggesting a high systemic inclination. 

We may have discovered a population of hidden accreting 
binaries. After optimally subtracting the stellar templates with 
matched spectral type, the residual spectra of seven sources show 
Ha emission. Three of them (CX138, CX377, and CXB73) have 
broad (> 400 km s“ ^) Ha lines, which are likely produced by an ac¬ 
cretion disk. CX84 and CXB82 are also possibly hidden accreting 
binaries based on their X-ray luminosity, while their residual Ha 
emission lines are narrow (<200 km s“^). Previous VLT/VIMOS 
spectra of CX377 showed a strong, double-peaked Ha emission 
line, while our Gemini/GMOS spectra of CX377 only contains an 
Ha absorption line. This may indicate the strong variability of the 
accretion disk of this source. The residual spectrum of CX377 after 
optimal subtraction supports this scenario by showing weak, broad 
Ha emission. 

In summary, based on the emission features, the spectral clas¬ 
sification, X-ray luminosity, and the residual spectra after optimal 
subtraction, we are able to constrain the likely nature for eight 
sources in our sample: CX446, CX1004, and CXB2 are accreting 
binaries (CVs or qLMXBs) while CX446 and CXB2 are likely also 
eclipsing; CXB64 and CXB113 are chromospherically active stars 
or binaries; CXI38, CX377, and CXB73 could be hidden accreting 
binaries. 


APPENDIX A: FINDING CHARTS 

The r'-band finding charts for the optical counterparts of 18 of the 
21 GBS sources presented in this paper are shown in Fig. [12] and 
Fig-El]. The sky area in each chart is 20" x 20" (except for CXB113 
which is 20" x 15") with the source in the center and indicated by 
the thick horizontal and vertical short bars. North is up and East is 
to the left. CXB 189 is blended with a nearby bright source because 
the field is crowded. Part of the emission from the nearby bright 
source may leak into the slit. The finding chart for CX377 is shown 
in Fig. [To] The charts for CX446 and CXI004 can be found in 
Torres et al. (2014). 
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Figure 12. The r'-band finding charts (generated from r'-band images obtained by the Blanco 4-meter telescope) for optical counterparts of 18 of the 21 
sources with Gemini/GMOS spectroscopy in this work (finding charts for the remaining three sources have been published in Torres et al. 2014; also see 
Fig.[l0|for the finding chart of CX377). The source is indicated by the short horizontal and vertical bars in each chart. The sky area of each chart is 20" x20", 
except for CXBl 13 which is 20" x 15". North is up and East is to the left. 
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Figure 13. Finding charts - Continued. 
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Table 1. The Observation Log of the Gemini/GMOS Spectroscopy 


GBS 

ID 

R.A. (J2000) 
(HHiMMiSS.SS) 

Dec. (J2000) 
(DD:MM:SS.S) 

Observation ID 
(GS-2012A-Q) 

Date (UT) 

Integration 
Time (s) 

Seeing 

(arcsec) 

CX84 

17:38: 12.84 

-29 : 06 : 12.4 

67-48 

2012-05-31 

1200 

0.5 




67-9 

2012-06-01 

1200 

0.5 

CX138 

17:46:23.14 

-25 : 49 : 30.3 

67-52 

2012-06-01 

1200 

0.5 




67-15 

2012-06-14 

1200 

1.1 

CX139 

17:45:22.15 

-25 : 50 : 48.2 

67-3 

2012-05-31 

900 

0.6 




67-44 

2012-06-01 

900 

0.6 

cx3n 

17 : 43 : 16.54 

-27 : 45 : 37.0 

67-70 

2012-06-01 

900x3 

0.6 




67-72 

2012-06-09 

900x3 

1.1 

CX446 

17:46:27.17 

-25 : 49 : 52.6 

67-96 

2012-06-22 

900x3 

0.8 




67-109 

2013-05-04 

900x3 

0.6 

CX1004 

17 : 46 : 23.47 

-31:05:49.8 

44-63 

2012-05-14 

900x3 

0.8 




44-59 

2012-05-17 

900x3 

0.8 




44-3 

2012-05-17 

900x3 

0.8 




44-61 

2012-05-17 

900x3 

0.7 

CXB2 

17 : 53 : 59.86 

-29 : 29 : 06.5 

44-53 

2012-04-20 

900x3 

0.9 




44-55 

2012-04-20 

900x3 

0.9 




44-57 

2012-04-20 

900x3 

0.8 




44-80 

2012-05-19 

900x3 

0.9 

CXB26 

17 : 53 : 47.89 

-29 : 44 : 37.4 

44-49 

2012-05-17 

900x3 

0.7 

CXB38 

17:58:33.81 

-27 : 30 : 22.9 

67-21 

2012-06-09 

900x3 

1.3 

CXB64 

17 : 46 : 02.87 

-32:08: 11.3 

67-74 

2012-05-29 

900x4 

0.7 

CXB73 

17 : 52 : 29.07 

-30:03:21.7 

67-68 

2012-06-22 

900x2 

1.1 

CXB82 

17 : 56 : 54.57 

-28 : 12 : 34.4 

67-66 

2012-06-14 

900x2 

1.0 




67-99 

2012-06-22 

900x2 

1.0 

CXB92 

17 : 34 : 25.02 

-30 : 41 : 04.6 

44-47 

2012-05-03 

900x3 

0.7 




44-85 

2012-05-19 

900 

1.1 

CXB99 

17 : 54 : 59.36 

-29: 10:21.1 

67-64 

2012-06-09 

900x2 

1.2 

CXB113 

17 : 52 : 05.53 

-30: 19:31.8 

67-62 

2012-06-22 

900 

1.1 

CXB117 

17 : 56 : 02.42 

-28 : 24 : 45.4 

67-60 

2012-06-01 

900x2 

0.6 

CXB137 

17 : 57 : 10.01 

-27 : 57 : 54.9 

44-45 

2012-05-17 

900x3 

0.8 

CXB149 

17 : 53 : 06.92 

-29:51 : 13.1 

67-58 

2012-05-31 

900x2 

0.5 

CXB174 

17:56: 18.52 

-28 : 45 : 40.5 

44-43 

2012-04-19 

900x3 

0.6 

CXB189 

17 : 52 : 55.98 

-29 : 29 : 50.7 

44-41 

2012-04-19 

900x3 

1.0 

CXB201 

17 : 33 : 26.08 

-30 : 40 : 24.4 

44-39 

2012-04-19 

900x3 

0.9 
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Table 2. Optical Photometry for the Gemini Counterparts 


GBSID 

Mean(r') 

Err(r') 

^ a(r') 

No. of Frames 

CX84 

19.018 

0.007 

0.029 

35 

CX377 

18.885 

0.009 

0.060 

36 

CX446 

21.162 

0.062 

0.166 

37 

CX1004 

20.761 

0.019 

0.040 

35 

CXB2 

20.350 

0.003 

0.14 

110 

CXB64 

- 18.63" 

0.011 

0.014 

56 

CXB82 

16.894 

0.006 

0.023 

52 

CXB99 

16.516 

0.021 

0.073 

110 

CXB113 

-19.13" 

0.018 

0.041 

55 


^The mean error of the r' magnitudes for the object. 

*^The standard deviation of the r' magnitudes for the object. 

^The magnitude here is approximate because of a lack of cali¬ 
brator stars in the field of view. 


Table 3. Near-Infrared Photometry for the Gemini Counterparts 


GBSID 

R.A.(VVV; J2000) 
(HHiMMiSS.SS) 

Dec.(VVV; J2000) 
(DD:MM:SS.S) 

Aa 

(arcsec) 

" J 

(mag) 

H 

(mag) 

K 

(mag) 

CX84 

17:38: 12.84 

-29 : 06 : 12.5 

0.12 

14.566 ±0.018 

13.570±0.019 

13.108±0.021 

CX138 

17:46:23.12 

-25 : 49 : 30.5 

0.36 

13.122 ±0.004 

12.134±0.004 

11.715 ±0.004 

CX139 

17:45:22.15 

-25 : 50 : 48.3 

0.22 

12.555 ±0.002 

11.914 ±0.003 

11.658 ±0.004 

CX377 

17 : 43 : 16.55 

-27 : 45 : 37.2 

0.20 

14.832 ±0.022 

13.870 ±0.026 

13.443 ±0.027 

CXB2 

17 : 53 : 59.85 

-29 : 29 : 06.4 

0.13 

17.251 ±0.344 

16.408 ±0.314 

16.251 ±0.336 

CXB38 

17:58:33.82 

-27 : 30 : 22.7 

0.25 

14.182±0.019 

13.264±0.016 

12.958±0.016 

CXB73 

17 : 52 : 29.08 

-30:03 :21.5 

0.18 

13.978 ±0.017 

13.143±0.016 

12.796 ±0.015 

CXB82 

17 : 56 : 54.57 

-28 : 12 : 34.3 

0.10 

13.514±0.011 

12.732 ±0.010 

12.536±0.011 

CXB92 

17 : 34 : 25.02 

-30 : 41 : 04.5 

0.10 

14.866 ±0.016 

13.967±0.018 

13.372±0.018 

CXB99 

17 : 54 : 59.34 

-29: 10:21.1 

0.19 

14.069 ±0.022 

13.356 ±0.022 

13.206 ±0.024 

CXB113 

17 : 52 : 05.52 

-30: 19:31.8 

0.10 

14.108 ±0.020 

13.599 ±0.023 

13.299 ±0.024 

CXB149 

17 : 53 : 06.93 

-29:51 : 13.0 

0.15 

14.242 ±0.022 

13.720 ±0.027 

13.527 ±0.029 

CXB174 

17:56: 18.53 

-28 : 45 : 40.5 

0.10 

14.388 ±0.029 

13.552 ±0.026 

13.268 ±0.025 


17 


The offset between near-infrared position and optical position, in arcseconds. 
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Table 4. List of Stellar Templates from the UVES POP Library 


Spectral Type Star Name 


AOV 

HD 162305 

AIV 

HD 65810 

A2V 

HD 60178 

A3V 

HD 211998 

A4V 

HD 145689 

A5V 

HD 39060 

A7V 

HD 187642 

A9V 

HD 26612 

FOV 

HD 109931 

FIV 

HD 40136 

F2V 

HD 33256 

F3V 

HD 18692 

F4V 

HD 37495 

F6V 

HD 16673 

F8V 

HD 45067 

F9V 

HD 10647 

GOV 

HD 105113 

GIV 

HD 20807 

G2V 

HD 14802 

G3V 

HD 211415 

G4V 

HD 59967 

G5V 

HD 59468 

G6V 

HD 140901 

G9V 

HD 25069 

K2V 

HD 22049 

K5V 

HD 10361 

MOV 

HD 156274 

M6V 

HD 34055 


Table 5. Speetral Classifieations for the Gemini Counterparts 


GBS ID Spectral Type 


CX84 

G9 

CX138 

G9 

CX139 

K2 

CX377 

F6 

CX1004^ 

M0-M5 

CXB64^ 

M0-M5 

CXB73 

G9 

CXB82 

G9 

CXB99 

K2 

CXB113^ 

M0-M5 

CXB149 

G6 

CXB174 

F4 


^These three sources can¬ 
not be spectrally classified via 
optimal subtraction. However, 
their spectra contain features 
(e.g., TiO band) of early-to- 
mid M-type stars. 
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Table 6. Radial Velocity Results for the Gemini Counterparts 


GBSID 

HJD 

CX84 

2456078.811 

2456079.677 

CX138 

2456079.760 

2456092.650 

CX139 

2456078.665 

2456079.783 

CX377 

2456079.714 

2456079.725 

2456079.736 

CXB73 

2456100.751 

2456100.762 

CXB82 

2456092.798 

2456092.809 

2456100.806 

2456100.817 

CXB99 

2456087.659 

2456087.670 


RVl (km RV2 (km s-^)^ 


0 

-104.3±3.3 

74.0±3.0 

-17.3±6.0 

0 

-110.7±2.1 

75.6 ±6.5 

-44.2 ±5.4 

0 

-96.2 ±10.6 

69.0±2.5 

-17.4±4.5 

0 

-20.5 ±9.2 

-14.6±12.9 

-43.6±11.9 

19.0±8.8 

-17.2±6.6 

0 

-30.0±6.2 

-1.7±5.1 

-31.0±6.0 

0 

-36.6±4.7 

-11.3±6.3 

-41.0±7.9 

-2.9±2.2 

-39.7±3.8 

-6.2±2.3 

-41.7±2.7 

0 

21.5±7.0 

-11.0±3.8 

-1.7±10.2 


Note. — Among the sources not listed in this table, seven of 
them (CX446, CX1004, CXB2, CXB64, CXB92, CXB113, and 
CXB137) do not show significant cross-correlations between their 
Gemini/GMOS spectra, while the other seven (CXB26, CXB38, 
CXB117, CXB149, CXB174, CXB189, and CXB201) have rela¬ 
tive radial velocities consistent with zero. 

^The radial velocity (RV) results when using the first source 
spectrum as a cross-correlation template; thus the RVl value for 
the first source spectra is always zero. 

’^The radial velocity (RV) results when using the star template 
spectrum with the spectral type shown in Table |5] as a cross¬ 
correlation template. 
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Table 7. Measurements of the Broad Ha Emission/Absorption Lines for Gemini Spectra 


GBS ID 

EW (A) 

FWHM(A) 

FWHM(kms-i) 

CX377 

3.1±0.6 

6.3 ±0.2 

233 ±9 

CX446 

-67.6 ±5.4 

27.5 ±1.2 

1250 ±50 

CX1004^ 

-38.0±0.6 

55.1 ±1.9 

2500 ±100 

CXB2 

-4.4 ±0.7 

18.6±3.7 

800 ±160 

CXB64 

-2.0±0.1 

4.9±0.1 

63±5 

CXB99 

-1.3±0.1 

6.3 ±0.6 

185 ±25 

CXB113 

-5.3±0.1 

5.3±0.1 

110±5 

CX84 

-2.8±0.1 

6.0±0.3 

225 ±15 

CX138 

-2.6±0.1 

30.0±1.5 

1350 ±75 

CX139 

-2.2±0.1 

6.4 ±0.6 

250 ±25 

CX377 

-2.8±0.1 

15.0±0.7 

660 ±30 

CXB73 

-1.3±0.1 

9.0±1.0 

350 ±50 

CXB82 

-1.3±0.1 

6.1 ±0.2 

180±10 

CXB99 

-4.1±0.1 

6.1 ±0.3 

180±15 


Note. — The upper part of the table shows the measurements of Ha 
emission/absorption lines in the Gemini/GMOS spectra of the GBS 
sources. The lower part of the tables lists the measurements of the 
Ha emission lines in the residual spectra. The residual spectra are 
obtained by optimally subtracting the best-ht stellar templates (see 
Fig-[I]for residual spectra). 

^The FWHM of the Ha emission/absorption lines before correcting 
for instrumental broadening, in unit of A. 

*^The FWHM of the Ha emission/absorption lines after correcting 
for instrumental broadening, in unit of km s“^. The FWHM values 
quoted in the text (in unit of km s“^) are also corrected for instrumen¬ 
tal broadening. 

^The measurements listed here are for the full Ha emission profile. 
For the measurements of each individual peak, see the text in ^4.1.11 








